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ABSTRACT 

We report the results of a high angular resolution near-infrared survey of dusty 
Wolf-Rayet stars using the Keck-1 Telescope, including new multi-wavelength im- 
ages of the pinwheel nebulae WR 98a, WR 104, and WR 112. Angular sizes were 
measured for an additional 8 dusty WR stars using aperture masking interfer- 
ometry, allowing us to probe characteristics sizes down to ~20 milliarcseconds 
(~40 AU for typical sources). With angular sizes and specific fluxes, we can 
directly measure the wavelength-dependent surface brightness and size relations 
for our sample. We discovered tight correlations of these properties within our 
sample which could not be explained by simple spherically-symmetric dust shells 
or even the more realistic "pinwheel nebula" (3-D) radiative transfer model, when 
using optical constants of Zubko. While the tightly-correlated surface brightness 
relations we uncovered offer compelling indirect evidence of a shared and dis- 
tinctive dust shell geometry amongst our sample, long-baseline interferometers 
should target the marginally-resolved objects in our sample in order to conclu- 
sively establish the presence or absence of the putative underyling colliding wind 
binaries thought to produce the dust shells around WC Wolf-Rayets. 
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1. Introduction 

The existence of dust shells around Wolf-Rayet stars has long posed a mystery: how can 
dust nucleate and survive near these hot stars in the presence of harsh ultraviolet radiation 
and gas densities lower than those found in other dust forming environments (e.g., AGB 
winds)? While novel chemical pathways for producing special carbonaceous grains have 
been proposed (Zubko 1998; Cherchneff et al. 2000), recent theoretical and observational 
results point strongly to a very different solution. 

Inspired by the periodic dust formation episodes around the colliding wind binary 
WR 140 (Moffat et al. 1987; Williams et al. 1990b), Usov (1991) suggested that dust for- 
mation might be catalyzed in the compressed layer of gas at the interface of colliding winds 
in WR+OB binary systems. This hypothesis was confirmed by the surprising imaging of a 
"pinwheel" nebula dust spiral around WR 104 (Tuthill et al. 1999b). Indeed, Monnier et al. 
(1999a) went so far as to suggest that perhaps all dusty WR systems are hiding colliding- 
wind binary systems (see also Dougherty & Williams 2000a) and this was further supported 
by the discovery of non-thermal radio emission around WR 104, WR 98a, and WR 112 
(Monnier et al. 2002b). 

The previously-cited imaging work was carried out using the diffraction-limited capabili- 
ties of the world's largest optical telescope at the Keck Observatory. High-resolution imaging 
was achieved using aperture masking interferometry, whereby the Keck-1 primary mirror is 
converted to a VLA-style interferometric array (Tuthill et al. 2000c) of many subapertures. 
This technique has been shown to be superior than current adaptive optics systems for 
high-fidelity imaging and size estimations for marginally-resolved objects (Rajagopal et al. 
2004). 

Here we report an extension to the initial imaging work of WR 104 and WR 98a. 
We observed all WR stars (accessible from Mauna Kea) brighter than m K 6.5 in order to 
investigate the binary hypothesis as to the origin of dust production. Although imaging 
was only possible for one additional source (WR 112), we report the characteristic size 
measurements for a total of 11 sources and discuss the significance of these findings. 

2. Observations 

Our group has been carrying out aperture masking interferometry at the Keck-1 tele- 
scope since 1996. We have published images and size measurements with unprecedented 
angular resolution on topics ranging from young stellar objects, carbon stars, red super- 
giants, and photospheric diameters of Mira variables (e.g., Monnier et al. 1999c; Tuthill 
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et al. 2000a,b; Danchi et al. 2001). A full description of this experiment can be found in 
Tuthill et al. (2000c), with further discussion of systematic errors in Monnier et al. (2004). 

The NIRC camera with image magnifier (Matthews et al. 1996) was used in conjunction 
with the aperture masking hardware. For this work, we used an aperture mask with an 
annulus that is 8-m in diameter as projected onto the Keck primary. This mask gives us 
complete UV coverage and sensitivity for targets brighter than itlk ~ 6 with a modest 
loss in calibration precision due to additional redundancy noise (for further discussion, see 
Tuthill et al. 2000c). The data frames were taken in speckle mode (T; nt =0.14s) to freeze 
the atmosphere. A variety of filters were used and the characteristic center wavelengths and 
widths can be found in Table 1 (spectral scans can be found in Keck-NIRC users manual). 

Here we report the full body of Wolf-Rayet size measurements collected during the 
entire period of the Keck aperture masking experiment. A target list and observing log can 
be found in Tables 2 & 3. All V 2 and closure phase data are available from the authors; all 
data products are stored in the FITS-based, optical interferometry data exchange format 
(OI-FITS), recently described in Pauls et al. (2005). 

3. Determining Characteristic Sizes 

3.1. Methodology 

While the basic data reduction has been described in previous papers (Tuthill et al. 
2000c), this paper is the first Keck aperture masking paper to explicitly deal with a sample 
of partially-resolved objects. This is important here since we intend to measure the 2- 
dimensional size and shape of Wolf-Rayet dust shells in order to search for asymmetries. If 
we measure an ellipticity in our data, what confidence do we have that an elongation is not 
due to poor calibration of the optical transfer function? To answer this question we have 
undertaken a systematic investigation of our calibration errors. 

The errors for this experiment are dominated by statistical and seeing calibration errors. 
The statistical error arises from the contribution of photon noise and read noise to our 
measurement process. Calibration error arises because the optical transfer function varies 
with time and telescope pointing. In order to correct for this latter effect, we always observe 
a point-source reference star nearby in time and angle from our target. In the process of 
calibration, raw visibilities from the target are divided by those from the calibrator. Hence, 
calibration error is multiplicative and affects high visibility data the most in absolute terms 
(AV). For well-resolved objects, these two types of error are comparable. For small objects 
(< 25 mas) calibration error is dominant and limits our ability to say whether an object is 
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resolved. 



3.1.1. Functional Form and Statistical Errors 

As part of our analysis, the calibrated visibility data are fit to the following generic 
function (based on a 2-dimensional Gaussian): 

V 2 = [V e -(- 2 W+7™) + v P }\ (1) 

where Vp is constant and represents a point source contribution (which was set to zero for 
this study). The parameters a, ft, and 7 are not constrained to be positive, thus allowing 
the function to curve upwards in cases where calibration errors are extreme. 

Under a rotation through angle 0, this function can be projected onto the major and 
minor axes (u',v'), written in the form 
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These expressions can be converted to the angular FWHM along the major/minor axis 
using the formula FWHM Major = and FWHM Minor = 

Statistical errors are dealt with by using bootstrap sampling (Efron & Tibshirani 1993) 
with the function given in Equation 1. The fitting routine MPFIT in Craig Markwardt's 
MINPACK 1 suite of IDL programs is used as the basis for the bootstrap. 



1 See http://cow.physics.wisc.edu/~craigm/idl/ 
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3.1.2. Monte-Carlo Method for Determining Calibration Error 

Quantifying the magnitude of calibration errors in this experiment has been difficult. 
Due to the very limited time available on a Keck telescope, the number of stars used as 
calibrators for a given filter is usually small, ranging from less than 5 up to perhaps a dozen 
for the most commonly used bands. While precise systematic errors are difficult to determine 
due to small number statistics, it is possible to estimate the magnitude of the calibration 
accuracy. 

The first step in the analysis of calibration errors is to cross-calibrate all of the calibration 
stars of a given night (with same bandpass filter and aperture mask) against each other and 
fit the resulting visibility curves with the quadratic function introduced in the last section 
(Eq.2). At this point, the calibrator-calibrator pairs are weighted according to proximity 

to each other in time and elevation according to e VAt * ; \ A6 ei*J with At* ~ 2000 s and 
A6 eU ~ 20°. 

In order to propagate calibration errors into the errors on the Gaussian FWHM sizes 
that we wish to measure, we start with a visibility curve with a known FWHM (using 
the same baselines as the Keck data). Then we multiply this perfect simulated data by 
various calibrator-calibrator pair data and then re-calculate a size using the quadratic model 
developed above. This simulates what the measured visibility would look like for a star of 
a known size using different calibrators. By Monte Carlo sampling over many cal-cal pairs 
according the weights described above, we build up an estimate of the likely error in the 
final size measurement. 

The results are summarized in Figure 1 for the 2000 June epoch in the CH4 filter using 
the Annulus mask. The most important feature of these plots is determination of the effective 
angular resolution limit of this experiment. We see that simulations of an unresolved point 
source, when passed through the calibration study described above, can yield apparent sizes 
as large as 12 milliarcseconds (Gaussian FWHM). From this, we conclude that our resolution 
limit corresponds to Gaussian FWHMs of 12-15 mas (A ~ 2.2/im), approximately 4x better 
than the formal Keck diffraction limit 1.22A/D ~ 55 mas. For nights of poor seeing, our size 
upper limit on FWHM increases to 15-20 mas. 

As stated before, we would like to use the Keck masking data to measure modest 
asymmetries in the dust shells of Wolf-Rayet stars, since this can tell us about the dust 
production mechanism (single-star vs binary interaction). In order to test our sensitivity 
to asymmetries, we fitted 2-D Gaussians to our calibrator study data and these results are 
also shown in Figure 1. Our study shows that if we impose typical calibration errors on a 
circularly symmetric object with FWHM 20 mas, we will regularly find asymmetries with 
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measured ratio FWHM minor / FWHM ma j or = 0.8. Thus, except for extremely asymmetric 
objects, any study on asymmetries from this data will have to restrict itself to objects that 
are larger than about >25-30 mas. 

We note that the calibration errors determined by this method are over-estimates, since 
typical calibrator-calibrator pairs suffer from greater time difference and angular distance 
than target-calibrator pairs used for actual science measurements (this is confirmed for a few 
objects where we have multiple independent data sets). Thus we can assume that the errors 
calculated in this section are conservative and the true calibration errors should be up to a 
factor of 1.5 to 2 smaller for nights of comparable seeing. 

The greatest uncertainty in this analysis is that it uses a very limited number of calibra- 
tor stars spread out over a night. Our understanding of the errors would improve greatly if 
a series of identical exposures of a calibrator star were to be taken sequentially 2 . This would 
allow us to see the short timescale variations in a calibrator's visibility curve, and thus how 
much variance there is in the calibration process. If very large variations in the quadratic fits 
remain on timescales as short as this, the calibration process may need further refinement. 

The same study has been done for H band (A = 1.65/zm) and PAHcs (A = 3.08/xm), and 
these results are also included (in the same format) in Figures 2 & 3. 

4. Results 

4.1. Multi- wavelength Characteristic Sizes 

Circularly-symmetric gaussian emission profiles were fitted to the azimuthally-averaged 
visibility data for all targets and the results can be found in Table 4. The uncertainty 
estimates were based on the combined error from statistical variations (using bootstrap) and 
from the calibration errors discussed in the last section. Calibration error is typically the 
dominant error for these data. 

The largest dust shells (WR 48a, WR 98a, WR 104, WR 112, and WR 140) show strong 
deviations from the Gaussian at the longer Keck baselines. This is not at all surprising since 
WR 98a, WR 104, WR 112 and WR 140 are all known to have spiral and/or asymmetric 
dust shells (Tuthill et al. 1999b; Monnier et al. 1999a; Ragland & Richichi 1999; Monnier 
et al. 2002d; Marchenko et al. 2002). It is likely that the other smaller dust shells in our 
sample are also not simple Gaussians, although we were unable to detect deviations due 



2 The Keck Time Allocation Committee is unlikely to give time to such a proposal 
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to insufficient angular resolution (all marginally resolved objects look Gaussian-like to an 
interferometer) . 

In order to calculate a consistent "characteristic size" for all targets, we only fit to 
short-baseline data where the visibility is greater than 0.5. By fitting only to short-baseline, 
high-visibility data, the reported characteristic sizes represent the overall scale of emission 
and are not strongly affected by the spiral and/or small-scale structures of the nebulae, if 
present. Our results are graphically shown for 2.2/im, 1.65//m and 3.08/zm in Figures 4-6. 

Subject to the (large) uncertainties in ellipticity discussed in §3.1.2, we did not find 
evidence for elongation in the emission except for the cases of the largest dust shells, WR 98a, 
WR 104, and WR 112. For these, we have reported the mean Gaussian FWHM and the ratio 
of the minor axis to the major axis. Just as for the characteristic sizes discussed above, this 
ellipticity parameter only refers to the "large-scale" emission component of the dust shell 
and not the smaller-scale details, such as the inner spiral windings of WR 104. 

4.2. Correction for Stellar Contribution 

While for most stars the near-infrared emission is completely dominated by dust emis- 
sion, some targets in our sample have significant stellar contributions. By carrying out simple 
spectral energy density (SED) fitting, we can estimate the fraction of light coming from the 
star (compared to dust) and apply a correction factor. Essentially the visible photometry 
is fit with a Kurucz model allowing the stellar size and reddening to vary; the stellar IR 
flux can then be estimated by extrapolation. With this knowledge, we can more accurately 
estimate the true size of the dust shell (not the star+dust emission). This is a common 
procedure in other areas of astronomy, such as estimating the sizes of disks around young 
stellar objects, and we follow procedures documented elsewhere (e.g., Millan-Gabet et al. 
2001; Monnier et al. 2005). 

This correction is strongest for 7 Vel (WR 11) which does not show much infrared 
excess, but since the observed characteristic size was consistent with a point source, we do 
not apply the correction. For cases where the correction is more than 5% (WR 48a, WR 95, 
WR 98a, WR 106, WR 113), we have included the new size estimates (and the estimates for 
the dust fraction) in Table 5. We estimate the correction factor is only known to 50% due 
to difficulty in uniquely fitting the broadband SEDs, and this error has been included in the 
corrected sizes in Table 5. 

Note that the post-outburst WR 140 dust shell is so large that it was obviously necessary 
to correct for the stellar contribution (see top curves of Figure 4) and this was done for 
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Table 4. The WR 140 dust shell was considered in detail in Monnier et al. (2002d) and we 
refer the reader to this paper for further discussion on this object. 

4.3. Aperture Synthesis Imaging 

Three systems in our sample are sufficiently resolved that they can be imaged using 
aperture synthesis techniques. The Maximum Entropy Method (MEM) (Skilling & Bryan 
1984; Gull & Skilling 1983) has been used to create diffraction-limited images from the 
interferometric data, as implemented in the VLBMEM package by Sivia (1987). 

Figures 7-9 show new multi- wavelength images of WR 98a, WR 104, and WR 112 
at wavelengths of 1.65/zm, 2.2/zm, and 3.08/im. The resolution degrades with increasing 
wavelength and we adopt an effective angular resolution of 21 mas, 28 mas, and 39 mas for 
1.65/zm, 2.2/zm, and 3.08/im, corresponding to A© = ^ — • While we have presented images 
of WR 98a and WR 104 in the past, images given here are at a new epoch and we have 
included additional wavelengths. 

Figure 9 contains the first published resolved images of the near-infrared dust shell 
around WR 112 (preliminary results were shown at the "Interacting Winds from Massive 
Stars" workshop in 2000; Monnier et al. 2002a). The asymmetric nature of this dust shell was 
already discovered using lunar occultations (Ragland & Richichi 1999). Recently, Marchenko 
et al. (2002) reported a spiral-like dust plume around WR 112 at mid-infrared wavelengths. 
Interestingly, we do not see an obvious spiral structure in these new near-infrared images, 
rather we see only evidence for filaments/arcs and one-sided nebulosity. 

5. Discussion 

This paper significantly expands the number of WR stars with known angular sizes, 
and we wish to use this information to probe the nature of the dusty outflows so common 
in WC systems. The Wolf-Rayet systems with the largest dust shells, WR 48a, WR 98a, 
WR 104, WR 112, WR 140, are all in confirmed colliding wind systems (from detection of 
non-thermal radio emission). Unfortunately, we do not have enough angular resolution to 
definitely resolve filamentary or spiral structure in the other objects in our survey. 
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5.1. Surface Brightness Relations 

Secure identification of binarity is critical - Monnier et al. (1999a) emphasize that if 
all dusty WR stars are in binaries that this has profound implications on massive stellar 
evolution and the nature of carbon-rich (WC) WR stars in particular. One method to 
achieve this end is to see if all the objects in our survey follow similar surface brightness 
relations. In order to calculate a specific surface brightness S\ we applied the following 
formula: 



F A lCH mag / 2 - 5 ) 

SX = O.68vr0| WHM (6) 

where F\ represent the flux density for a zero magnitude star (here using units of W 
m~ 2 /im _1 ), mag is the magnitude of the WR dust shell at the given wavelength (correcting 
for stellar emission when appropriate - see §4.2), and Ofwhm is the FWHM of the Gaussian 
fit. Note the 0.68 in the formula is a correction factor which converts a Gaussian FWHM 
size estimate to the equivalent uniform disk for a better definition of "Surface Brightness." 
The final units of Sa are W m" 2 /im _1 sr _1 . 

Figure 10 shows the observed surface brightness relations. The objects in the survey 
all appear to follow the same surface brightness relations, indicating all the systems share 
a similar near-infrared emission mechanism. Assuming the emission comes from optically- 
thin carbon dust (described by the optical constants of Zubko 1998), we can estimate the 
color temperature, and we find T co i or ~ 1000 K (considering the relation between 1.65/im 
and 2.2/zm) and T co i or ~ 650 K (considering 2.2/xm and 3.08/im). The fact that the surface 
brightness relations can not be described by a single temperature reflects the fact that the 
dust from a wide range of temperatures contributes to the near-IR emission. 

We compared our observed surface brightness relations with the predicted relations from 
two models. Firstly, we used the radiative transfer code DUSTY (Ivezic et al. 1999) to fit 
spherically-symmetric dust shell models (assuming uniform outflow and dust constants of 
Zubko 1998) to the SED. We could vary the amount of dust optical depth and calculate 
the model surface brightness subject to an overall fit to the SED. We present the calculated 
surface brightness relations in Figure 10, spanning T2.2^ =0.01-0.77. We note that the SED 
fits were not perfect (usually under-predicting 3.1/zm fluxes), and thus we take these relations 
as representative but not fully optimized. 

Harries et al. (2004) recently calculated the emergent flux from a series of spiral dust 
shell models using the 3-D (Monte Carlo) radiative transfer code TORUS (Harries 2000). We 
refer to this as the "WR 104 Reference Model" and used this sophisticated physical model to 
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calculate the surface brightness relations as a function of inclination angle (the model only 
reasonably fits the WR 104 SED for low inclinations). These results can also be found on 
Figure 10. 

Comparing the model calculations to the data, we find that neither of these models 
are a good fit to all the 1.65/im, 2.2/iin, and 3.1//m surface brightness relations, although 
the models can reproduce some of the observations. More generically, this figure shows us 
that any particular dust shell geometry yields a family of surface brightness relations. The 
tightly-correlated surface brightness relations measured for our sample is suggestive of a 
shared and distinctive dust shell geometry, although neither the spherically-symmetric nor 
3-D WR 104 Reference Model can fully explain our results. One could improve the fits of a 
given model to the data by tuning the adopted optical constants (which are not well known), 
although the differences between models would persist. Before discussing the implications of 
this, we look for more information in the "size ratio" relations. 



5.2. Size Ratio Relations 

Another possibly distinctive (and distance-independent) signature to differentiate spherically- 
symmetric dust shells from colliding wind systems is the ratio of dust shell characteristic sizes 
at different wavelengths. We have collected all of the size ratios and presented these (with 
errors) in Figure 11. Assuming all the objects are drawn from the same class, we find the 
following mean size relations: 

FWHM at 1.65/xm 

^1.65 M m/2.2 M m = ™ mA/r . „ = 0.73 ± 0.05 (7) 



FWHM at 2. 2 urn 

FWHM at 3.08/xm 
iW»/2.*™ - pwHM at 2 2fxm ~ 1-36 ±0.07 (8) 



This compares favorably with calculations for an optically-thin, spherically-symmetric 
dust shell (-Ri.65/im/2.2 M m = 0.78, -R 3 .08Mm/2.2 Aim = 1-22), as performed using DUSTY (Ivezic 
et al. 1999) and tuned to fit the SED. We have also calculated the size ratios predicted by the 
WR 104 reference model discussed in the last section. As shown in Harries et al. (2004), this 
model does a remarkable job in fitting the morphology, SED, and overall size. The series of 
dust shell models used to calculate the size relations (as well as surface brightness relations) 
can be found in Figure 12. Analysis of these synthetic images yields the following mean size 
relations: 



inclination 0° 



-Rl.65/um/2.2/im 



= 0.92 



- 11 - 



inclination 60' 



inclination 30' 



inclination 90' 



R3.08tim/2.2fj,m = 1-04 

Rl.65nm/2.2nm = 0.91 

R3.08fim/2.2fj,m = 1-07 

Rl.65fim/2.2fmi = 0.95 

R3.08fim/2.2f J ,m = 1-14 

R\m\xtnj2.2iim = 0.98 

R3.08fmi/2.2fj,m = 1-28 



Paradoxically, the simple spherically-symmetric model fits the WR 104 size relations 
better than the model made specifically for WR 104. This can be explained if the optical 
depth of the shell in the reference model of Harries et al. (2004) is somewhat too high. 
High r in the inner windings of the spiral casts a strong shadow across the nebula which 
cause the temperature profile to be very steep (analogous to the temperature profile for an 
optically-thick dust shell). Most likely the observed size relations can be accommodated 
by the pinwheel model by adjusting the distribution of dust in the inner most part of the 
nebulae. 



WR 112: We know that WR 112 is in a colliding wind binary system (like WR 98a 
and WR 104) from its clear non-thermal radio emission (Chapman et al. 1999; Monnier et al. 
2002b). Furthermore, we believe that we are viewing the binary near the orbital plane (not 
face-on) because the non-thermal radio emission is highly variable, suggesting that our line- 
of-sight periodically passes through the optically-thin O-star wind (as is well-documented 
for WR 140, White & Becker 1995). In this context, it is not surprising that the colliding 
wind dust spiral takes on a more complicated shape due to projection effects (see discussion 
of this in Monnier et al. 2002d; Harries et al. 2004). Indeed, the % =60° model for WR 104 
(Figure 12) bears an obvious similarity to the "horseshoe" structure that we observed for 
WR 112 (middle-panel of Figure 9). 

This suggestion of an edge-one viewing angle contradicts the apparent (near) face-o dust 
spiral geometry seen in the mid-IR (Marchenko et al. 2002). We speculate that projection 
effects of the outer spiral windings (seen at high inclination) could also explain the mid- 
IR data, yielding bright arcs and filaments that may resemble a face-on spiral outflow (e.g., 
Monnier et al. 2002b). Alternatively, if the underlying colliding wind system is indeed viewed 
from a face-on position as implied by the mid-IR images, then the variable radio emission 



5.3. Comments on individual objects 
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might be caused by a, as yet undetected, third stellar component to the system. Further 
analysis of this system will be included in a future paper where we will present multi-epoch 
near- and mid-IR images as well as VLA monitoring photometry. 

WR 113 CV Ser: van der Hucht (2001) lists WR 113 as a nearly edge-on binary 
system (inclination~70°) with colliding winds. We expect such edge-on systems to show 
significant dust shell elongations in near-IR (see Figure 12) , although our data indicates the 
dust emission is symmetric (Ratio r = FWHM minor /FWHM ma j or >0.8). This higher level of 
symmetry could result if the cone opening angle of colliding wind interface was larger than 
found for WR 104 or if there was significant dust entrained in the WR outflow in a more 
spherical pattern. Perhaps the hint of a 3rd component in the system (Niemela et al. 1996) 
could also explain this surprising lack of dust shell asymmetry in this object. 

WR 118: Our K-band size estimate is compatible with the earlier speckle measurement 
of Yudin et al. (2001). 

Quintuplet Stars: Tuthill et al. (2006) recently reported resolving two bright Quin- 
tuplet WC stars into pinwheel nebulae. These objects show a size ratio Rzm^m/2.2nm ~ 2, 
larger than for the stars reported here, along with a correspondingly low 3.1/im surface 
brightness. We have not included these stars in this paper given the difficulty in correcting 
for extinction and the possible importance of local heating in this unusually dense and active 
cluster. 



6. Conclusions and Future Work 

We have presented a multi-wavelength survey of near-infrared angular sizes of dusty 
Wolf-Rayet systems from the Keck Aperture Masking program. Aperture synthesis images 
were presented for WR 112 for the first time and we found strong evidence for interact- 
ing/colliding winds. In addition, we presented new epochs of WR 104 and WR 98a images 
as well as the first results at 3.1/im, confirming the spiral nature of the dust shells for WR 104 
and WR 98. 

Using these data, we discovered tightly-correlated surface brightness relations and also 
common size ratios between different near-infrared bands. The observed relations could not 
be reproduced in detail using either a spherically-symmetric dust shell model or the 3-D 
radiative transfer model of WR 104 (Harries et al. 2004). The high-quality data presented 
here will act as an observational foundation for a new generation of modelling efforts. 

We find these results to be compelling indirect evidence that all these dusty WR here 
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share a common emission geometry, presumably related to the obvious spiral distribution of 
WR 104 and WR 98a (and confirmed recently for two Quintuplet WC stars; Tuthill et al. 
2006). Thus, while not conclusive, this study can be viewed as further evidence that WC 
Wolf-Rayet stars are associated with binarity (Monnier et al. 1999a; Williams & van der 
Hucht 2000). 

In order to make further progress on identifying binarity, we must pursue multiple 
approaches. While traditional direct spectroscopic identifications of binarity have proved 
difficult (e.g., Williams & van der Hucht 2000), other promising methods have emerged for 
unambiguously determining if dusty WR systems are binaries. We recommend the following 
strategies be pursued: 

• Higher resolution infrared data. A factor of 3 longer baselines will allow unambiguous 
detection of asymmetry and non-zero closure phases, if spirals and/or filaments exist 
in our sample. This corresponds to about 30 m baselines, easily achievable with the 
current interferometers (the Very Large Telescope is particularly well-suited to this). 
The main challenge will be the low visible flux which makes tip-tilt tracking difficult - 
infrared star trackers would be helpful in this regard. 

• Visible or infrared photometry. Photometric fluctuations have been reported for WR 98a 
(Monnier et al. 1999a) and WR 104 (Kato et al. 2002). The variations seem to correlate 
to the known orbital periods and we strongly encourage observers to monitor all dusty 
WC stars to establish the periods of the putative underlying binaries. This method 
can work for objects at large distances (e.g., galactic center, Local Group), however 
one has to contend with intrinsic variability and have the determination to monitor 
sources over many years (or decades; Williams et al. 1990a). 

• Non-thermal radio emission. For systems with periods ^1 year, we expect detectable 
non-thermal radio emission from the colliding winds (Dougherty & Williams 2000b; 
Monnier et al. 2002c). This requires sensitive measurements at multiple frequencies 
but should be possible and will yield unambiguous results. Recent successful attempts 
to use this method include Leitherer et al. (1997), Chapman et al. (1999), Monnier 
et al. (2002b), and Cappa et al. (2004). 

While skeptics will not be convinced yet that the WC phenomena has a necessary con- 
nection to binarity, the evidence in favor of this scenario continues to accumulate. The 
proposed observational efforts to conclusively establish binarity complement on-going theo- 
retical studies to explore the far-reaching consequences of the role of binarity in the theory of 
massive star evolution. Indeed, a complete understanding of Gamma Ray Burst progenitors, 
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the Galactic Black Hole population, and the high-mass end of the Initial Mass Function 
hinges on an accurate picture of massive binary stellar evolution. 
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Table 1. Properties of NIRC Camera Infrared Filters 



Name 


Center Wavelength 


Bandpass FWHM 


Fractional 




A (H 


AA (/im) 


Bandwidth 


Fell 


1.6471 


0.0176 


1.1% 


H 


1.6575 


0.333 


20% 


K 


2.2135 


0.427 


19% 


Kcont 


2.25965 


0.0531 


2.3% 


CH4 


2.269 


0.155 


6.8% 


PAHcs 


3.0825 


0.1007 


3.3% 



Note. - Reference: The NIRC Manual (Harrison & Goodrich 
1999) 



Table 2. Basic Properties of Targets 



Source 


RA (J2000) 


Dec (,12000) 


V 


J 


H 


K 


Spectral 


Distance 


Lum 


inosity 


Photometry 


Names 














Type 


(kpc) 


Log I 


:L/L©) b 


References 


7 2 Vcl, WR 11 


08 09 31.96 


-47 20 11.8 


1.81 


2.15 


2.25 


2.10 


WC8 + 07.5III-V (1) 


0.26 (2) 


5.5° 


(3,4) 


5 


WR 48a 


13 12 39.65 


-62 42 55.8 




8.74 


6.80 


5.09 


WC8ed + ■>. (1) 


1.2 (1) 






5,9 


WR 76 


16 40 05.3 


-45 41 10 


15.36 


8.46 


6.51 


4.88 


WC9d (1) 


5.35 (1) 






5,9 


WR 95 


17 36 19.76 


-33 26 10.9 


14.00 


8.29 


6.67 


5.27 


WC9d (1) 


2.09 (1) 






5,9 


WR 98a 


17 41 13.051 


-30 32 30.34 




9.14 


6.51 


4.33 


WC8-9vd + OB? (1) 


1.9 (8) 






9 


WR 104 


18 02 04.123 


-23 37 42.24 


13.54 


6.67 


4.34 


2.42 


WC9d + B0.5V (1) 


2.3 (6) 


5.4° 


(7) 


5 


WR 106 


18 04 43.66 


-21 09 30.7 


11.93 


7.94 


6.28 


4.82 


WC9d (1) 


2.3 (1) 






5,9 


WR 112 


18 16 33.489 


-18 58 42.47 


17.7 


8.68 


6.26 


4.26 


WC9d +OB? (1) 


4.15 (1) 






5,9 


WR 113 


18 19 07.36 


-11 37 59.2 


9.43 


7.02 


6.28 


5.49 


WC8d +08-9 (1) 


1.79 (1) 






5,9 


WR 118 


18 31 42.3 


-09 59 15 


22 


8.10 


5.41 


3.65 


WC9 (1) 


3.13 (1) 






5,9,10 


WR 140 


20 20 27.98 


+43 51 16.3 


6.9 


5.55 


5.43 


5.04 


WC7pd + 04-5 (1) 


1.85 (11) 


6.1° 


(12) 





a These magnitudes (V band from Simbad, and J.H.K bands from 2MASS) are merely representative since some of the targets are variable. 

b Luminosity here has been corrected for the adopted distance used in this paper; not all objects have well-established luminosities due to high extinction. 

c Binary system total luminosity. 



Note. — References: (1) van dcr Hucht (2001), (2) ESA (1997), (3) Do Marco et al. (2000), (4) do Marco & Schmutz (1999), (5) Williams et al. (1987), (6) , 
Tuthill et al. (1999a), (7) Harries et al. (2004), (8) Monnicr et al. (1999b), (9) 2MASS; Cutri et al. (2003), (10) Simbad, (11) Dougherty et al. (2005), (12) 

Williams et al. (1990a) ^ 

I 
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Table 3. Journal of Observations 



Target 


Date 


Filter (A ) 


Aperture 


Calibrator 41 




(UT) 


(Mm) 


Mask 


Names 


7 2 Vol 


1999 Feb 05 


CH4 (2.269) 


Golay-21 


HD 68553 (K3, 4.39±0.44 mas) 


w r!, 4tta 


zUUU J an zo 


IS. [Z.ZlOO ) 


Annulus-36 


Mil llooyy ^rio, l.z±0. / masj 


WE 76 


1998 Jim 05 


H (1.6575) 
K (2.2135) 


Annulus-36 


HD 151834 (G3, 1.2±0.4 mas) 




1999 Apr 26 


H (1.6575) 
K (2.2135) 

.FAriCS {o.VQZo) 




HD 153258 (K4, 2.4±0.9 mas) 


WE 95 


2000 Jun 24 


H (1.6575) 
K (2.2135) 


Annulus-36 


HD 158774 (M0, 2±2 mas) 


WE 98a 


2000 Jun 24 


H (1.6575) 

~LT fO Ol QK 1 ) 

IS. \Z . zloo ) 
CH4 (2.269) 
PAHcs (3.0825) 


Annulus-36 


HD 158774 

Hd 163428 (K5, 4.1±1.0 mas) 


WE 104 


2000 Jun 24 


H (1.6575) 
CH4 (2.269) 
PAHcs (3.0825) 


Annulus-36 


HD 165813 (M0, 1.7±0.9 mas) 
HD 167036 (K2, 3.0±0.6 mas) 
HD 167036 


WE 106 


1998 Jun 05 

1999 Jul 30 


H (1.6575) 
CH4 (2.269) 
PAHcs (3.0825) 
K (2.2135) 
PAHcs (3.0825) 


Annulus-36 


HD 164124 (K2, 1.1±0.3 mas) 
HD 167036 
HD 167036 


WE 112 


2000 Jun 24 


H (1.6575) 
K (2.2135) 
CH4 (2.269) 
PAHcs (3.0825) 


Annulus-36 


HD 165813 
HD 167036 


WE 113 


2000 Jun 24 


H (1.6575) 
K (2.2135) 


Annulus-36 


HD 168366 (K2, 1.03±0.12 mas) 
HD 168366 


WE 118 


1998 Jun 05 


H (1.6575) 
CH4 (2.269) 
PAHcs (3.0825) 


Annulus-36 


HD 168000 (K3, 1.42±0.8 mas) 
HD 175775 (K0, 3.5±1.1 mas) 




1998 Jun 06 


H (1.6575) 
CH4 (2.269) 
PAHcs (3.0825) 


Annulus-36 


HD 170474 (K0, 1.3±0.3 mas) 
HD 173074 (M, 3.4±1.4 mas) 
HD 175775 




1999 Apr 26 


H (1.6575) 
CH4 (2.269) 
PAHcs (3.0825) 


Annulus-36 


HD 168000 
HD 175775 


WE 140 


1999 Jul 30 


Fell (1.6471) 
Kcont (2.25965) 


FFA 
FFA 


HD 193631 (K0, 0.5±0.1 mas) 




2001 Jul 30 


Fell (1.6471) 
Kcont (2.25965) 


FFA 
FFA 


HD 192867 (Ml, 1.6±0.6 mas) 






PAHcs (3.0825) 


FFA 


HD 192867 



a Calibrator size estimates made using getCal, maintained and distributed by the Michclson Sci- 
ence Center (http://msc.caltech.edu). 
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Table 4. Mean Infrared Size Results 



Tarj 


ret 


Date 


Gaussian FWHM ( 


mas) a,i> 


Comments 






(UT) 


1.65/iin 


2.2/im 


3.08/um 




7 2 Vel 


1999 Feb 05 




7±7 






WR 


48a 


2000 Jan 26 




63±3( c > 




Fit to "large-scale" component only. 


WR 


76 


1998 Jun 05 


17±3 


19±4 










1999 Apr 26 


19±3 


20±4 


21±4 




WR 


95 


2000 Jun 24 


16±3( c > 


23±4( c > 






WR 


98a 


2000 Jun 24 


73±4( c > 


110±7 


147±9 


Fit to "large-scale" component only. 








r=0.56 


r=0.74 


r=0.81 


Elliptical (PA~120°). 


WR 


104 


2000 Jun 24 


57±9 


73±7 


101±9 


Fit to "large-scale" component only. 








r=0.84 


r=0.88 


r=0.85 


Elliptical (PA~30°). 


WR 


106 


1998 Jun 05 


17±3( c > 


23±4( c > 


33±3 








1999 Jul 30 




21±4( c > 


24±4 




WR 


112 


2000 Jun 24 


40±10 


83±7 


131±9 


Fit to "large-scale" component only. 










r=0.87 




Elliptical (PA~130). 


WR 


113 


2000 Jun 24 


19±2( c > 


27±3( c > 






WR 


118 


1998 Jun 05 


18±3 


24±4 


30±3 








1998 Jun 06 


18±3 


23±4 


32±3 








1999 Apr 26 


21±3 


23±4 


33±3 




WR 


140 


1999 Jul 30 


<15 


<13 




Pre-outburst 






2001 Jul 30 


268±17 


228±12 


263±13 


Two-component model: Gaussian + Point 








39% dust 


64% dust 


78% dust 





a These characteristic sizes are not corrected for the presence of the star (which contribute significant near-IR 
flux in some cases). 

b If visibility data are manifestly elliptical, then we report both the mean FWHM and the Ratio r = 
FWHM minor /FWHM major . 

c See Table 5 for corrected sizes. 



-22 - 



Table 5. Characteristic sizes of Wolf-Rayet dust shells, corrected for stellar contribution 



Target Dust Fraction (%) Gaussian FWHM (mas) a 
1.65/im 2.2/im 1.65/im 2.2/xm 



WR 48a 




85% 




70±5 


WR 95 


73% 


89% 


19±4 


25±4 


WR 98a 


90% 




78±5 




WR 106 


78% 


93% 


20±4 


24±5 






93% 




22±5 


WR 113 


38% 


65% 


35±10 


35±6 



a These characteristic sizes are corrected for the pres- 
ence of the star (which contribute significant near-IR flux 
in some cases; see §4.2). Refer to Table 4 for size measure- 
ments for rest of sample which do not require corrections. 



- 23 - 



i 





¥ 




r 

¥ 











10 20 30 40 50 60 

True FWHM (mas> 




0.2 - 

o.o r, 1 1 1 1 1 1 1 1 1 1 1 

10 20 30 40 50 

True FWHM (mas> 



Fig. 1. — These graphs show the results for the fitting of simulated data using calibrators 
from the June 2000 epoch for the CH4 band using the Annulus mask. The top panel shows 
the measured Gaussian FHWM from the simulated data compared to the input model. 
The middle panel re-displays this information as a ratio. The bottom panel explores how 
miscalibration induces ellipticity into the fits, even for an input circular Gaussian. The 
dashed lines indicate the position of the original simulated data. All error bars indicate 1-a 
certainty and data points indicate medians. 



True FWHM (mas) 



True FWHM (nasi 



Fig. 2. — These graphs show the results for the fitting of recalibrated simulated data using 
calibrators from the June 2000 epoch for the H band using the Annulus mask. The format 
of the plots are the same as previous figure. 
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Fig. 3. — These graphs show the results for the fitting of recalibrated simulated data using 
calibrators from the June 2000 epoch for the PAHcs band using the Annulus mask. The 
format of the plots are the same as previous figure. 
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Fig. 4. — This figure shows the azimuthally-averaged visibility curves at 2.2yum for Wolf- 
Rayet stars in our sample, along with best-fit Gaussian curves. Each visibility curve is 
offset by 0.25 from the next. For the most resolved dust shells, the Gaussian curves were 
only fitted to the "large-scale" component of the visibility data (V > 0.5). See Table 4 for 
detailed fitting results. 
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Fig. 5. — This figure shows the azimuthally-averaged visibility curves at 1.65/im for Wolf- 
Rayet stars in our sample, along with best-fit Gaussian curves. Each visibility curve is 
offset by 0.25 from the next. For the most resolved dust shells, the Gaussian curves were 
only fitted to the "large-scale" component of the visibility data (V > 0.5). See Table 4 for 
detailed fitting results. 
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Fig. 6. — This figure shows the azimuthally-averaged visibility curves at 3.1/xm for Wolf- 
Rayet stars in our sample, along with best-fit Gaussian curves. Each visibility curve is 
offset by 0.25 from the next. For the most resolved dust shells, the Gaussian curves were 
only fitted to the "large-scale" component of the visibility data (V > 0.5). See Table 4 for 
detailed fitting results. 



WR 98a 
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Contours (% of Peak): 1 2 4 8 16 32 64 

Fig. 7. — Here we show multiwavelength aperture synthesis images of WR 98a on UT 2000 
Jun 24, at 1.65//m, 2.2/im, and 3.08/im. Note that the resolution is lower at the longer 
wavelengths, as indicated by the "beam" spot located in the bottom-left corner of each 
panel (representing the best achievable angular resolution, A© = — )• 




Fig. 8. — Here we show multiwavelength aperture synthesis images of WR 104 on UT 2000 
Jun 24, at 1.65/xm, 2.2/im, and 3.08/im. Note that the resolution is lower at the longer 
wavelengths, as indicated by the "beam" spot located in the bottom-left corner of each 
panel (representing the best achievable angular resolution, A© = — )• 
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Fig. 9. — Here we show multiwavelength aperture synthesis images of WR 112 on UT 2000 
Jun 24, at 1.65/xm, 2.2/im, and 3.08/im. Note that the resolution is lower at the longer 
wavelengths, as indicated by the "beam" spot located in the bottom-left corner of each 
panel (representing the best achievable angular resolution, A© = ^ — )• 
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Fig. 10. — Here we show the surface brightness relations for dusty Wolf-Rayet systems. We 
plot the best-fit linear relation corresponding to T co i or ~ 1000 K (left panel) and T co i or ~ 
650 K (right panel), assuming optical constants for Zubko (1998) dust. The open circles 
show the surface brightness relations for a series of spherically-symmetric dust shell models 
with r 2 . 2At =0.01, 0.1, 0.3, 0.67 (higher r yields higher S\). The filled circles show the same 
relations for the WR 104 Reference Model of Harries et al. (2004) for inclination angles of 
% =0, 30, 60, 90° (higher % yields higher S\). The arrow shows how the surface brightness 
changes for interstellar reddening of Ay = 5, assuming standard ISM dust. 
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Fig. 11. — Here we plot near-infrared characteristic size ratios vs K-band magnitude. The 
left panel shows the results for a * t , while the right panel shows the results for 

p^^m a ^ t 3 2 ° 2 8 ^ . We plot (dashed line) the mean size ratio and also mark with arrows the ex- 
pected size ratios for (a) optically-thin, spherically-symmetric dust shell and (b) the WR 104 
"Reference Model" of Harries et al. (2004) viewed from 30° inclination angle. 
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Fig. 12. — Near-IR Size Study of WR 104 Reference Model. This figure shows a gallery of 
synthetic images of the WR 104 reference model described in Harries et al. (2004), smoothed 
to a uniform angular resolution of 20 mas. These images were used to derive the surface 
brightness and size ratio comparisons described in §5.1 & 5.2. The 3 columns correspond to 
effective wavelengths of 1.65/im, 2.2/im, and 3.08/xm, while the 4 rows correspond to viewing 
angles of inclination 0°, 30°, 60°, and 90°. 



